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Abstract 

We discuss the possibility of accurately estimating the source number density of 
ultra-high-energy cosmic rays (UHECRs) using small-scale anisotropy in their ar- 
rival distribution. The arrival distribution has information on their source and source 
distribution. We calculate the propagation of UHE protons in a structured extra- 
galactic magnetic field (EGMF) and simulate their arrival distribution at the Earth 
using our previously developed method. The source number density that can best 
reproduce observational results by Akeno Giant Air Shower Array is estimated at 
about 1CT 5 Mpc -3 in a simple source model. Despite having large uncertainties of 
about one order of magnitude, due to small number of observed events in current 
status, we find that more detection of UHECRs in the Auger era can sufficiently 
decrease this so that the source number density can be more robustly estimated. 200 
event observation above 4 x 10 19 eV in a hemisphere can discriminate between 10~ 5 
and 10 -6 Mpc -3 . Number of events to discriminate between 10 -4 and 10~ 5 Mpc -3 
is dependent on EGMF strength. We also discuss the same in another source model 
in this paper. 
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1 Introduction 

The nature of the sources of ultra-high-energy cosmic rays (UHECRs) is still 
poorly known though many efforts in detection of these particles with UHE 
energies. It is one of challenging problems in astroparticle physics. 

1 E-mail addresses: takami@utap.phys.s.u-tokyo. ac.jp (H. Takami), sato@phys.s.u- 
tokyo.ac.jp (K.Sato) 



Preprint submitted to Elsevier Science 



2 February 2008 



About 10 years ago, Akeno Giant Air Shower Array (AGASA) reported small- 
scale anisotropy of the UHECR arrival distribution within its angular resolu- 
tion while it did large-scale isotropy with a harmonic analysis [1,2]. It found 
five doublets and one triplet as event clusterings. It is enough evidence that 
origin of UHECRs is astrophysically point-like source. However, no obvious 
astronomical counterparts to the observed UHECRs have been found. One of 
this reasons is magnetic fields in the universe. 

Extragalactic magnetic field (EGMF) is poorly known theoretically and ob- 
servationally. As an upper limit for its strength, B, and correlation length, 
l c , B l^ 2 < (1 nG)(l Mpc) 1 / 2 , by Faraday rotation measurements of radio 
signals from distant quasars, is often adopted [3]. Based on its upper limit, 
UHECR propagation has been discussed in simply uniform turbulent magnetic 
field [4,5,6]. 

It is also observationally known that clusters of galaxies have strong magnetic 
fields with 0.1 - a few /iG at its center [7]. This is shown that there are not 
only uniform EGMF but also relatively strong EGMF incidental to the local 
structure. Recently, several groups have performed simulations of large-scale 
structure formation with magnetic fields [8,9,10]. They find that magnetic field 
traces the local density distribution. Such magnetic field plays an important 
role in UHECR propagation since it is stronger than the uniform one. 

In our previous study [11], we constructed a structured EGMF model that 
reflects the local structures and discussed predicted arrival distribution of UHE 
protons at the Earth in consideration with their propagation processes. We also 
considered Galactic magnetic field (GMF) in the propagation. As a result, the 
number density of UHECR sources was constrained to ~ 10~ 4 Mpc -3 under 
our luminosity-weighted source model (explained in section 2) from observed 
arrival distribution. The EGMF strength was normalized only to 0.4 fxG at the 
center of the Virgo cluster. However, observational measurements of magnetic 
field in a galaxy cluster result in 0.1 - a few fiG, so that the strength of 
magnetic field has large uncertainty of about one order of magnitude. Thus, 
it is important that the propagation process and the arrival distribution are 
investigated in several strengths of EGMF. 

In this study, we calculate propagation of UHE protons in a structured EGMF 
with several strengths for normalization and simulate the arrival distributions 
at the Earth. At first, comparing them to observational results, we constrain 
the source number density. We use AGASA data for this purpose because 
AGASA has observed the most number of events. Such constraint is dependent 
on UHECR source model. We discuss about two simple source models. We find 
that the number density has large uncertainty of about an order of magnitude 
due to small number of observed events at present. So, we also discuss the 
possibility of a decrease in the uncertainty with future observations. 
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In section 2, our models of UHECR source distribution, structured EGMF, and 
GMF are briefly explained. In section 3, a numerical method of UHE proton 
propagation, construction of the arrival distribution and statistical analysis 
are explained. We present our results in section 4, and summarize in section 
5. 



2 Our Models 

In this section, our distribution models of UHECR sources and magnetic fields 
are briefly explained. These models are almost the same as those in our pre- 
vious work, thus, more detailed explanation is written in reference [11]. 

Our models of the source distribution and a structured EGMF are constructed 
out of the Infrared Astronomical Satellite Point Source Catalogue Redshift 
Survey (IRAS PSCz) catalog of galaxies [12]. This catalog is thought to be 
the most appropriate galaxy catalog for these purposes since it has the largest 
sky coverage (~ 84% of all the sky). The selection effects in the observation are 
corrected with a luminosity function of the IRAS galaxies [13]. A set of galaxies 
after this correction is called our sample galaxies below. Our sample galaxies 
within 100 Mpc are used for construction of our models since small number 
of galaxies can be detected outside 100 Mpc. We adopt fl m = 0.3, Q\ = 0.7 
and H = 71 km s _1 Mpc -1 as the cosmological parameters. 

We assume a subset of our sample galaxies to be an UHECR source distribu- 
tion. The number density of UHECR sources is taken as our model parameter. 
For a given source number density, we randomly select galaxies from our sam- 
ple galaxies with probabilities proportional to the absolute luminosity of each 
galaxy. Source distribution outside 100 Mpc is assumed to be isotropic and 
luminosity distribution of the galaxies follows the luminosity function. 

In this study, we adopt two simple source models. One is a source distribution 
that all sources have the same power for injection of UHE protons, called nor- 
mal source model, and the other is that each source has a power proportional 
to its luminosity, called luminosity-weighted source model. In both the models, 
it is assumed that all sources have the same of maximum acceleration energy. 
Our sample galaxies are over about five orders of magnitude in the luminosity. 
We construct 100 source distributions in each source number density. 

Several simulations of the large-scale structure formation with magnetic field 
have found that the EGMF roughly traces distribution of the baryon density 
[8,10]. According to these results, our structured EGMF model is constructed 
based on simple assumptions. We constructed the matter density distribution 
from our sample galaxies with a spatial resolution of IMpc. In each cell, a 
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strength of the EGMF is related to the matter density, p, as |-E>egmf| oc p 2 / 3 
and a turbulent magnetic field with the Kolmogorov spectrum is assumed. 
The strength of EGMF is normalized at a cell that contains the center of the 
Virgo Cluster to 0.0,0.1,0.4 and 1.0/iG. Note that magnetized space covers 
about 5% of the universe within 100 Mpc. The remaining 95% space has no 
magnetic field. Outside of 100 Mpc, we assume an uniform turbulent field with 
1 nG. 

As GMF model, a model in reference [14] is adopted. This model consists of a 
bisymmetric spiral field and a dipole field. We neglect a turbulent component 
of the GMF though its field strength is in the range of 0.5-B rcg — 2B vcg where 
B TCg is the strength of the regular component of the GMF. The turbulent 
component little changes the arrival direction of UHE protons [15]. 



3 Numerical Methods 



3. 1 Calculation of the Arrival Distribution of UHE Protons 



In this study, propagation of UHE protons is calculated by an application 
for the backtracking method. UHE protons injected from their sources cannot 
always reach the Earth because of their deflections by magnetic fields. It wastes 
many CPU time to calculate the propagation of such protons and to construct 
the arrival distribution. In order to solve this problem, we suggested a new 
method of calculation of the propagation in intergalactic space with attention 
to an inverse process of UHECR propagation [11]. In this method, we consider 
UHECRs with the charge of -1 that ejected from the Earth, and calculate their 
trajectories in Galactic space and intergalactic space. Trajectories of these 
particles can be regarded as trajectories of UHE protons from extragalactic 
space. 

In intergalactic space, we consider not only the deflections due to EGMF, 
but also the energy-loss processes, which are adiabatic energy loss due to the 
expanding universe, the electron-positron pair creation and photopion produc- 
tion in collision with the Cosmic Microwave Background (CMB) [16,17]. All 
of these processes are treated as continuous processes. 

The adiabatic energy loss results from the expansion of the universe. This 
energy loss rate is written as 



dE a 

— = — E = -H Q 

at a 



1/2 



n m (i + z y + n A e. (i) 
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For the pair creation we adopt the analytical fitting functions given by [18] to 
calculate the energy-loss rate on isotropic photons. For the photopion produc- 
tion, the energy-loss length which is calculated by simulating the photopion 
production with the event generator SOPHIA [19] is used. 

In Galactic space, all energy loss processes are neglected since the path lengths 
of UHE protons are enough shorter than energy-loss lengths of those processes. 

The arrival distribution of UHE protons at the Earth can be generated from 
their trajectories calculated with a method explained above. We eject 2,500,000 
particles from the Earth isotropically with the ejection spectrum of dN/dE oc 
E~ l '° and record their trajectories. The trajectories are calculated until their 
propagation times reach the age of the universe or their energies reach 10 22 eV, 
which is assumed as maximum energy of accelerated proton at each source. 
With our source models, we calculate a factor for each trajectory, which rep- 
resents the relative probability that the j th proton reaches the Earth, 



P rFoWV" hhi dN/dEg(dij,E g l ) dE g 

Pselec(^j) OC ^ — — — . (2) 



Here % labels sources on each trajectory, while Zij, dij, and L^j are the red- 
shift, distance, and luminosity of such sources. E g is the energy of a proton 
at a source, which has the energy E at the Earth. Thus, dN / dE g (d i ^ ) E g l ) is 
the energy spectrum of UHE protons ejected from a source whose distance is 
dij. We assume that UHE protons are injected with a power-law spectrum in 
the energy range of 10 19 — 10 22 eV. dE g /dE represents a correction factor for 
the variation of the shape of the energy spectrum through the propagation. 
The power-law index of dN/dE g is assumed to be 2.6 in order to fit the cal- 
culated energy spectrum above 10 19 eV to the observational result [20]. This 
corresponds to an injection model that the ankle is assumed to be the pair 
creation dip [21]. 

We randomly select trajectories according to these relative probabilities, -P S eiec, 
so that the number of the selected trajectories is equal to the required event 
number. The mapping of the ejected direction of each particle at the Earth is 
the arrival distribution of UHE protons. If we have to select the same trajec- 
tory more than once to adjust the number of selected trajectories, we generate 
a new event whose arrival direction is calculated by adding a normally dis- 
tributed deviation with zero mean and variance equal to the experimental 
resolution to the original arrival direction. 
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3.2 Statistical Methods 



In this subsection, we introduce two statistical quantities needed to compare 
the arrival distribution of UHE protons with that obtained by the observation 
statistically. 

The two-point correlation function N(9) contains information on the small- 
scale anisotropy. We start from a set of events generated from our simulation. 
For each event, we divide the sphere into concentric bins of angular size A9 
and count the number of events falling into each bin. We then divide it by the 
solid angle of the corresponding bin, that is, 



where denotes the separation angle of the two events. The angle A9 is taken 
to be 1° in this analysis. The AGASA data show a correlation at small angles 
(~ 3°) with a 2.3(4.6)<x significance of deviation from an isotropic distribution 
for E > 10 19 eV(E > 4 x 10 19 eV) [2]. 

In order to investigate what number density of sources can best reproduce the 
two-point correlation observed by AGASA, we introduce Xe m ^ f° r a source 
distribution as 



Xe max 



^ [N(0) - N^(0)Y ^ (4) 



where N(9) is the two-point correlation function calculated from a simulated 
arrival distribution within —10° < 5 < 80°, which corresponds to the sky 
observed by AGASA and N ohs (9) is that obtained from the AGASA data at 
angle 9. The o~(9) is the statistical error of N(9) due to the finite number of 
simulated events. The random event selection is performed 20 times. This X6» max 
represents the goodness of fit between the simulated two-point correlation and 
the observed one. In this study, we take 9 max to be 10°. 
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Fig. 1. xios as a function of the source number density in the normal source model. 
The error bars originate from 100 times source selection. The GMF is considered in 
the lower panel while not in the upper panel. 

4 Results 

4-1 Constraints on the Source Number Density 

The arrival distributions of UHE protons are constructed by a method ex- 
plained in section 3.1. Comparing the simulated arrival distribution with the 
observed arrival distribution statistically, we constrain the number density of 
UHECR sources. 

First, we try to restrict the number density from small-scale anisotropy. As 
mentioned above, the two-point correlation function is a statistical indicator 
of the small-scale anisotropy. Therefore, we compare the two-point correlation 
functions calculated from simulated arrival distributions and the functions 
calculated by AGASA result above 4 x 10 19 eV [22]. In recent year, the au- 
thors of reference [23] pointed out that small-scale anisotropy by AGASA is 
statistically not consistent with the observed spectrum above 10 20 eV (above 
Greisen-Zatsepin-Kuz'min(GZK) cutoff [24,25]) by them. Our source models 
also predict GZK cutoff even if our structured EGMF and GMF are taken into 
account, as shown in next subsection. Thus, we assume that the super-GZK 
events are of another origin. So, in this study we restrain the energy range of 
observed events we use up to 10 20 eV. 

Figure 1 shows Xio between two-point correlation functions calculated in the 
normal source model and AGASA data as a function of the number density of 
UHECR sources. The upper panel shows Xio s calculated without the GMF for 
B = O.OyuG (crosses), B = 0.1/iG (triangles), B = 0.4/xG (squares), and B = 
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Fig. 2. The two-point correlation functions calculated from only source distributions 
that can predict the large-scale isotropy. The histograms are the observational result 
within 4 x 10 19 < E < 10 20 eV (49 events) by AGASA. The error bars are from the 
event selection for finite events and the shaded regions show total la errors. The 
GMF is included. 



1.0/xG (pentagons) respectively. The error bars come from cosmic variance due 
to 100 source distributions for every number density. The lower panel shows the 
same, but GMF is included. The same number densities are represented despite 
four marks being a little shifted to the horizontal axis for easy visualization 
to see. 



The smaller xw gives the better number density. Xios of 10~ 7 Mpc~ 3 are 
relatively larger than those of the others. Hence, 10~ 7 Mpc~ 3 cannot best 
reproduce the AGASA result. However, the remaining x 10 s provide similar 
values. More discussion is needed in order to estimate the number density 
more precisely. 
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Fig. 3. Same as fig. 1, but calculated in the luminosity-weighted source model. 

UHECR source distribution must satisfy not only the small-scale anisotropy 
but also the large-scale isotropy. The large-scale isotropy is estimated quan- 
titatively by harmonic amplitudes. Thus, we recalculate the two-point cor- 
relation function only from source distributions that predict the large-scale 
isotropy comparable to AGASA, which is quantified using the same method 
as in reference [11]. 

Figure 2 shows such two-point correlation functions with several values of the 
number density and strengths of the EGMF. The GMF is also considered. The 
histograms are the two-point correlation functions calculated from AGASA 
events with energies from 4 x 10 19 to 10 20 eV (49 events). The error bars 
are due to the finite number of event selection, and the shaded regions are 
total lo" error including errors due to different source distributions. The errors 
that originate from different source distributions are very small since source 
distributions that predicts large values of the function are excluded. Such 
source distributions cannot satisfy the large-scale isotropy. 

We focus on the left three panels with Beg — 0.0/iG. In the upper panel, 
the mean values (points) are smaller than the observational data and cannot 
reproduce it within la error at small angular scale. On the other hand, source 
distributions with number densities of 10~ 6 Mpc -3 predict stronger anisotropy 
at the small angle scale even if the large-scale isotropy is satisfied. Thus, 
the appropriate number density that best reproduces the arrival distribution 
observed by AGASA is ~ 1(T 5 Mpc" 3 in the case of .Beg = 0.0/iG. This is 
the same result as reference [26,27] though they do not consider the GMF. 

In the other panels, similar tendencies to the case of B^g — 0.0/iG can be 
found. The EGMF diffuses cosmic rays during propagation, thereby weakening 
event clusterings. Therefore, smaller number densities (~ 10~ 6 Mpc -3 ) are 
favored in stronger EGMFs. However, in the case of 1(T 6 Mpc -3 with B EG = 
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Fig. 4. Same as fig. 2, but calculated in the luminosity-weighted source model. 

0.4, 1.0/iG, number of source distributions used is only about 30% of 100 source 
distributions while about 70% at 10~ 5 Mpc~ 3 . Although a stronger EGMF 
allows a smaller number density, the appropriate number density is close to 
10~ 5 Mpc -3 . 

From the discussion above, the number density of 10~ 5 to 10~ 6 Mpc~ 3 can 
best reproduce the AG AS A results, but, as you find in the upper panels, with 
uncertainty of about one order of magnitude. The result is almost unchanged 
even if the GMF is neglected. 

We can discuss similarly in the case of the luminosity weighted source model. 
Figure 3 shows the same as figure 1, but calculated with the luminosity- 
weighted source model. The error bars are larger than those in figure 1 because 
of additional degree of freedom in the UHECR source luminosities. Mean val- 
ues of Xw & re smallest at 10~ 4 Mpc~ 3 , about one or two orders of magnitude 
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larger than those in the case of the normal source model, but with some un- 
certainty. Thus, we discuss the arrival distributions that satisfy the large-scale 
isotropy again. Figure 4 shows the same as figure 2, but calculated based on 
the luminosity-weighted source model. Two-point correlation functions pre- 
dicted from the source number density of ~ 1CT 3 Mpc~ 3 are smaller than 
the observational results at small scales. In the lower left panel, the predicted 
function is too large to reproduce the observed one at 2 — 4° within la er- 
rors. When the strengths of the EGMF are stronger, values of the functions 
at small angles are more suppressed and the two-point correlation functions 
reproduce the observational result well. However, the number of source distri- 
butions with 10 -5 Mpc -3 , adopted in this calculation, is only about 10% while 
the number with 10" 4 Mpc~ 3 is about 50%. Hence, the appropriate number 
density is thought to be closer to 10 -4 Mpc -3 . 

As a result, in the luminosity-weighted source model, the number density of 
UHECR sources is 10~ 4 Mpc~ 3 with uncertainty. This is about one order of 
magnitude larger than that in the normal source model. As mentioned in 
section 2, the injection luminosity of sources is over 5 orders of magnitude in 
the luminosity-weighted source distribution. In these sources, more luminous 
sources can more strongly contribute to the arriving cosmic rays. Therefore, 
the number of such contributed sources seems smaller than the real number 
of cosmic ray sources. 



4-2 Energy Spectra above 10 19 eV 

We discuss the dependency of the energy spectra on the magnetic field and 
source models. Figure 5 shows the energy spectra calculated from 100 source 
distributions with 10 -5 Mpc -3 in the normal source model (left panel), and 
10 -4 Mpc -3 in the luminosity-weighted source model (right panel) . These spec- 
tra are predicted in the cases of B-^q = 0.0 (black), 0.1 (red), 0.4 (green), and 1.0 
/j,G (blue). The error bars are la total statistical error including cosmic vari- 
ance. The GMF is also considered. The spectra are normalized at 10 19 - 55 eV. 
The histogram show the AGASA results [28] with shaded region representing 
statistical errors. 

Around 10 19 eV, a stronger EGMF predicts a little lower flux since it lengthens 
the propagation distance and UHE protons lose their energies more readily. 
Under this normalization, a relatively higher flux is predicted at the highest 
energy region in a stronger EGMF. However, these are not statistically signif- 
icant. Our structured EGMF cannot affect the energy spectrum. Our source 
models also cannot reproduce the super-GZK events. We assume that the 
super-GZK events are of other origins. 
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Fig. 5. Energy spectra of UHE protons calculated from all of 100 source distribution 
with 10 -5 Mpc -3 in the normal source model (left panel), and from that with 
10~ 4 Mpc~ 3 in the luminosity- weighted source model (right panel). Color points with 
error bars, which originates from la total statistical error, are calculated spectra 
with 0.0 (black), 0.1 (red), 0.4 (green), and 1.0 /t/G (blue). They are normalized at 
10 19 ' 55 eV. The histogram and shaded region are AGASA results [28]. 



4-3 Future Prospects of our Approach 



In section 4.1, the number density of UHECR sources that can best reproduces 
the arrival distribution observed by AGASA is constrained. However, it has 
a large uncertainty that originates from the small number of observed events. 
Thus, our next interest is whether increase in observed events can determine 
the number density more accurately or not. In order to discuss this problem, 
it is necessary to compare our predicted arrival distribution to the arrival dis- 
tribution observed by future observations. However, because we cannot know 
future results, we adopt an isotropic distribution of UHECR arrival directions 
as a template for future observations. A similar approach was adopted in ref- 
erence [29]. However, simulation of [29] did not include magnetic fields. We 
discuss the future prospects taking into account GMF and EGMF. 

The small-scale anisotropy will be more strongly observed if UHECRs come 
from some point-like sources. Therefore, we compare the two-point correlation 
functions of our simulated arrival distribution to that calculated from the 
isotropic arrival distributions. To make the comparison, we introduce \ 2 as 

N bi n e=1 o CT sim (9) 2 + CT iso (9) 2 ' 



where N sim (9) and N iso (9) are the two-point correlation functions calculated 
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Fig. 6. Distributions of y^s, calculated from arrival protons above 4 x 10 19 eV, 
at several strengths of the EGMF in the normal source model. The GMF is not 
considered. The strengths of the EGMF are 0.0fiG(upper left), 0.lfiG(upper right), 
OAfiG(lower left), and 1.0/xG (lower right). The numbers of events are set to be 49 
events within —10° < 5 < 80°, 200 events and 500 events within the southern 
hemisphere to emulate Auger. 

from simulated and isotropic arrival distributions respectively, and <r sim (9) and 
Ciso(^) are la statistical errors of N sim (9) and N iso (9) due to the finite number 
of events. iVbin is set to be 5 to obtain information on the small-scale anisotropy. 
X 2 s are distributed with their averages and variances for every source number 
density. If these distributions can be distinguished, we can argue that the 
source number density can be determined more accurately. 

Figure 6 shows distributions of x 2 in the normal source model at several EGMF 
strengths. The strengths of the EGMF are 0.0, 0.1, 0.4, and 1.0// G. These are 
calculated from arrival protons above 4 x 10 19 eV. The numbers of events are 
set to 49 events within -10° < 5 < 80° to emulate AGASA, and 200 events 
and 500 events within —90° < 5 < 0° to simulate Pierre Auger Observatory 
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(Auger) [30]. The GMF is not considered, but the GMF hardly affects the 
result. The separation of the \ 2 distribution from the first bin represents the 
appearance of difference from isotropic event distribution. 

At current status (49 events), the \ 2 distribution is localized at x 2 — 1 f° r 
1(T 4 and 1(T 5 Mpc" 3 . This shows that 70-80% of the source distributions 
with such number densities have an isotropy consistent with the large-scale 
isotropy reported by AGASA. The distributions with 10~ 6 Mpc~ 3 is a little 
shifted to the right for _B EG = 0.0/iG. This shows that the arrival distribution 
expected for 10~ 5 Mpc -3 can be distinguished from that of 10~ 6 Mpc~ 3 , to 
some extent. This is also supported in fig. 1. The others with 10~ 6 Mpc~ 3 
overlap the distribution for larger densities. Thus, the distinction between the 
number densities is difficult at present. It is remarkable that a stronger EGMF 
predicts a distribution biased to the left since it has the arrival distribution 
to be more isotropic. 

Next, we discuss a near future status with 200 events above 4 x 10 19 eV, 
which corresponds to the number of events observed by Auger up to 2007. 
The distributions for 10 -6 Mpc~ 3 do not overlap with those of the other num- 
ber densities in all panels. Therefore, Auger can observe a difference between 
10" 6 Mpc~ 3 and more number densities while it is difficult to discriminate 
10~ 5 Mpc -3 from 10 -4 Mpc~ 3 except for the case of Beg — 0.0/xG. Source 
distributions with less than 10~ 5 Mpc -3 do not predict isotropy anymore. 

Accumulation of 500 events gives rise to larger separation. Even the arrival 
distributions expected from 10 -4 Mpc~ 3 are distinguishable from isotropic 
distribution. For Beg — 0.0/xG, the three distributions can be completely sep- 
arate, and, we can determine the source number density from an observational 
arrival distribution. However, when the EGMF is considered, the distributions 
with 10" 5 Mpc~ 3 overlap with those of 10" 4 Mpc -3 , since the stronger EGMF 
diffuses the UHECR arrival directions and their arrival distribution becomes 
more isotropic. More observation is required. 

In the luminosity-weighted source model, we can discuss as above for 10~ 3 , 10~ 4 
and 10~ 5 Mpc. The graph is shown in figure 7. At present (49 events), 10~ 5 Mpc 
has distributions with large width. This fact reflects that only 10% of the 
source distributions can reproduce the large-scale isotropy. For 200 event ob- 
servation, the three distributions are separated more largely. In this case, one 
question occurs. The source distributions with 10~ 3 Mpc -3 in Beg — 0.1/xG 
do not predict an isotropy similar to that in Beg — 0.0/xG. This seems to 
be strange. This is because the arriving cosmic-rays from near sources rel- 
atively increase at highest energies by the EGMF. Source distributions with 
10~ 3 Mpc -3 contain the most number of UHECR sources. Because near sources 
cause anisotropy, isotropy is worse than that in Beg = 0.0/xG. When the 
EGMF is stronger, the diffusion of cosmic-rays can overcome the anisotropy. 
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Fig. 7. Same as fig.6, but the luminosity-weighted source model is adopted. 



Hence, the distributions of x 2 f° r B EG = 0.4 and 1.0/iG are shifted to the left. 
200 event observation enables us to distinguish 10~ 3 Mpc" 3 with 10~ 4 Mpc -3 . 
The arrival distribution from less than 10~ 4 Mpc -3 is distinguishable from the 
isotropic distribution. 



For 500 event observation, 10~ 3 Mpc" 3 is separated from the others because 
of strong deflection providing UHECR arrival distributions isotropy. However, 
the distributions of 10" 4 Mpc -3 and 10" 5 Mpc" 3 have broad width due to the 
degree of freedom of the source luminosity. Even 500 events observation cannot 
enable us to discriminate these source number densities, but to discriminate 
10" 3 Mpc" 3 from the isotropic distribution. 
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5 Summary & Discussion 



In this study, we discuss the possibility of accurately estimating the source 
number density of UHECRs with small-scale anisotropy. Comparison between 
simulated arrival distribution and the observational results enables us to esti- 
mate the source number density. In order to construct the arrival distribution, 
we calculate the propagation of UHE protons in a structured EGMF with 
several strengths consistent with measurements of magnetic field in clusters 
of galaxies. The GMF is also considered. We find that the source number 
density of 1CT 5 Mpc~ 3 in the normal source model and 10~ 4 Mpc -3 in the 
luminosity-weighted source model can best reproduce the AGASA results, 
which are weakly dependent on strength of our structured EGMF. However, 
these have large uncertainty of about one order of magnitude due to the small 
number of observed events. 

So, we discuss the possibility that future observations decrease the uncer- 
tainty. In the normal source model, Auger can distinguish 10~ 5 Mpc~ 3 and 
1(T 6 Mpc~ 3 sufficiently by our method in the near future. If the structured 
EGMF is zero or very weak, 10~ 4 Mpc -3 is also discriminated from the 
less number density in 500 event observation above 4 x 10 19 eV. In stronger 
EGMF, more observations are requested because cosmic rays are diffused more 
strongly. Number of events that needed for the distinction depends on EGMF 
strength. In the luminosity-weighted source model, 10 -3 Mpc~ 3 can be dis- 
tinguished from the less number density by Auger. The distinction between 
the less number densities is difficult due to large uncertainty which originates 
from different injection powers of the sources. 

We exclusively adopt the AGASA results in this study. However, High Reso- 
lution Fly's Eye (HiRes) claims no significant small-scale anisotropy, contrary 
to AGASA [31]. This discrepancy is one of well-known problems in UHECR 
experiments. At present, this is not statistically significant due to the small 
number of observed events [32]. It will be able to solved by new experiments 
with large aperture, such as Auger, Telescope Array [33], and Extreme Uni- 
verse Space Observatory [34]. Hence, it is possible that these experiments 
do not observe sufficient small-scale anisotropy in the future. If Auger does 
not observe small-scale clusterings during 2007 (maybe it detects about 200 
events above 4 x 10 19 eV), the source number density is estimated at about 
10 -4 Mpc~ 3 or more in the normal source model by definition of x 2 i n figure 
6. It is comparable with number density of active galactic nuclei[35]. 

Our EGMF model within 100 Mpc has about 95% of volume without magnetic 
field. In our model, uniform magnetic field is not considered. According to an 
upper limit mentioned in section 1, deflection angle of UHE protons with 
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Fig. 8. Distribution of x 2 s calculated from arrival protons above 4 x 10 19 eV, in the 
normal source model. The GMF is not considered The strength of the structured 
EGMF is normalized to 0.1 fiG and that of an uniform turbulent field is (upper 
panel), 1 nG (middle panel), and 10 nG (lower panel). The numbers of events are 
set to be 200 (left) and 500 events (right) within the southern hemisphere. 

energy of E during propagation of distance, d, is estimated as 
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This deflection is expected to generate more isotropic arrival distribution of 
UHECRs. Therefore, uniform turbulent magnetic field affects determination 
of the source number density. 

As a demonstration, we show x 2 distribution for B EG = 0.1 /iG, including the 
uniform field, in figure 8. The normal source model is adopted. The number 
of events is set to be 200 (left panel), and 500 (right panel). The strengths 
of the uniform turbulent field are (upper panel), 1 (middle panel), and 10 
nG (lower panel). The distributions are shifted to lower value in stronger 
uniform turbulent field. The shifts are larger in more source number density. 
In stronger turbulent uniform magnetic field, number of events needed for 
the distinction between 10~ 5 and 10~ 6 Mpc -3 is smaller while it becomes 
difficult to discriminate 10~ 4 from 10~ 5 Mpc~ 3 due to the diffusion of cosmic 
rays. Strength of uniform EGMF is also important for estimating the source 
number density. 

In this work, we adopt isotropic arrival distribution of UHECRs as a template 
of a future result since results of new experiments are still unpublished. Auger 
should detect a few times more number of UHE events than that of AGASA 
since it started observation about three years ago. Its result will provide us 
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beneficial information on the nature of UHECR sources. 
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